Relations between the diffusion anomaly and cooperative rearranging regions in a 
hydrophobically nanoconfined water monolayer 
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We simulate liquid water between hydrophobic walls separated by 0.5 nm, to study how the 
diffusion constant Dv, parallel to the walls depends on the microscopic structure of water. At low 
temperature T, water diffusion can be associated with the number of defects in the hydrogen bond 
network. However, the number of defects solely does not account for the peculiar diffusion of 
water, with maxima and minima along isotherms. Here, we calculate a relation that quantitatively 
reproduces the behavior of Du , focusing on the high-T regime. We clarify how the interplay between 
breaking of hydrogen bonds and cooperative rearranging regions of 1 nm size gives rise to the 
diffusion extrema in nanoconfined water. 

PACS numbers: 68.08.-p, 66.10.C- 68.15.+e 



In normal liquids the diffusion constant D decreases 
when the pressure P increases at constant temperature 
T. Water, instead, displays up to a 60% increase of D 
[l[ for P < 200 MPa . Based on simulations of bulk 
water with detailed models, e.g. Refs. and lattice 

models, e.g. relations between the behavior of D and 
the structure of water have been proposed. For example, 
D can be related to the configurational entropy Q, and 
the minima in D can be associated with a maximum in 
orientational order @. In classical molecular dynamics 
(MD) simulations an increase of P weakens the hydrogen 
bonds (HBs), thus increasing D This interpretation 
in terms of defects in the HB network can be extended 
to negative P [5j. A similar qualitative conclusion has 
been reached also by ab initio MD Q. Nevertheless, a 
quantitative relation between the microscopic structure 
of water and the anomalous behavior of D, including both 
extrema, is still missing. 

Diffusion experiments in hydrophobic confinement are 
at issue. For water confined in carbon nanotubes (CNTs) 
with diameters from 5 nm to 2 nm, D decreases for 
smaller diameters [10|. Nevertheless, other experiments 
reveal an exceptionally fast transport for water confined 
in CNTs of about 2 nm [ll[ and 7 nm diameter fl2j . 

Results from models are also controversial. MD sim- 
ulations of five point transferable intermolecular poten- 
tial (TIP5P) water nanoconfined between hydrophobic 
smooth walls display anomalous diffusion constant D\\ 
parallel to the walls at lower T than in bulk [l3j], and 
no anomaly in the orthogonal direction [TiJ . A large de- 
crease of diffusion is found at ambient conditions for ex- 
tended simple point charge (SPC/E) water between two 
large hydrophobic graphite-like plates for separations be- 
low 1.3 nm [la ]. Nevertheless, first-principles MD of the 
same model in similar conditions show that the diffu- 
sion becomes faster under confinement, possibly due to 
weaker HBs at the interface (lrlj . A similar controversy 
is reported for simulations of water in CNTs, with diam- 
eters below 1 nm 



Here we study, by the Monte Carlo (MC) method, the 
diffusion of water in a monolayer between two parallel hy- 
drophobic walls, adopting an already well-studied coarse- 
grained model [19| . We find maxima and minima of Dm 
and derive an expression that relates D\\ to water config- 
urations, identifying cooperative rearranging regions and 
a relation among structure, thermodynamics and dynam- 
ics that clarifies the mechanisms for water diffusion. 

We consider two flat hydrophobic walls at a distance of 
h ~ 0.5 nm such that the formation of ice is inhibited |20j. 
Adopting the natural square symmetry of the system [20j , 
we divide the available volume V into J/ square cells, 
each with a volume v — V j ,JV , and hydrate the system 
with N < JV water molecules. To each cell we associate 
an occupation variable = 0,1 (i = 1,2,..., ,A') if it 
is vacant or occupied, respectively. The enthalpy of the 
system is 



PV, 



(M)« 



where ry is the distance between water molecules i and 
j, U(r) = oo for r < tq = 2.9 A, the water van der 
Waals diameter, U(r) = e w [(r /r) 12 — (r /r) 6 ] for r > r 
with e w = 5.8 kJ/mol, the isotropic attraction energy, 
and U(r) = for r > r c = \f~Jf /4, the cut-off dis- 
tance. J = 2.9 kJ/mol is the characteristic energy of 
the directional (covalent) component of the HB |2l| . 
A^hb = E(i j) n i n j^dij,(jji is the total number of HBs, 
where the sum is over nearest-neighbor (n.n.) molecules, 



1 if o-j 



otherwise, and the 



variables cry are defined as follows. We adopt a geomet- 
rical definition in which the HB breaks if OOH > 30°. 
Therefore, only 1/6 of the orientation range [0,360°] 
in the OH-0 plane is associated with a bonded state. 
Hence, we account for the entropy loss of molecule i due 
to the formation of a HB with molecule j by associat- 
ing to i a bonding index cry, and to j a bonding index 
<7ji, with both indices € [1,2, ... ,q] with q = 6. For the 



square symmetry, each molecule has four n.n. and four 
dij, with q 4 = 6 4 = 1296 possible bonding states. There- 
fore, when two n.n. molecules i and j with mrij = 1 form 
a HB, their energy and entropy decrease. 

HB formation increases the average volume per 
molecule, because it induces a local structure with a re- 
duced number of n.n. as compared to close packing [22j . 
This effect is accounted for by an enthalpy increase Pvhb 
for each HB, where uhb/^o = 0.5 is the average density 
increase from low density ice Ih to high density ices VI 
and VIII, and v$ = hr^. The total volume occupied by 
water is V w = Nv + -/Vhb^hb ■ The increase «hb cor- 
responds to a larger volume per molecule, but not to a 
larger separation r between molecules, having no effect 
on U(r). 

In Eq. (JTJ), Jo- = 0.29 kJ/mol is the energy gain for two 
bonding indices of the same molecule in the same state, 
and accounts for the HB many-body (cooperative) inter- 
action 0, HH, with the sum over (k,l)i, i.e., over each 
of the six different pairs of the four er^ of a molecule i. 
This cooperative interaction among HBs favors specific 
T and P dependent values of the probability distribu- 
tion of 0-0-0 angles 23]. In confinement, the distri- 
bution displays no fifth interstitial n.n. and has a max- 
imum shifted toward 90° at low T [24J, consistent with 
the square symmetry adopted here. The water-wall in- 
teraction is represented by a hard-core exclusion. 

We perform MC simulations in the NPT ensemble 
for ,JV — 2500 and N/,yY — 0.75%, corresponding to 
N = 1875 water molecules. Since we allow for changes of 
the volume in the direction parallel to the walls, the con- 
trol parameter P represents the pressure parallel to the 
walls. We test that for jY — 400 and 1600 at the same 
75% occupancy ratio there are no appreciable differences, 
as well as for 75% < N '/ ' JY < 90%. The detailed MC 
algorithm and the conversion to real units are described 
elsewhere [n|. We equilibrate each state point for 0.2 ms 
and average over 15 ms [25]. 

We calculate Dm from the two-dimensional Einstein re- 
lation 

D„ = km <l*(' + '°) -*(*)!») 

11 t-VOG At ^ ' 

where r 4 (t) is the projection onto the plates of the posi- 
tion of molecule i at time t, and (•) is the average over 
all molecules i and over different values of to. We find 
that D|| decreases below a maximum P/j" ax for decreasing 
pmax < q 2 GPa, as in bulk water Q, reaching a mini- 
mum D™ in at P min < P max (Fig. H) . P min decreases with 
decreasing T and eventually becomes negative. For the 
high temperatures considered here, D\\ is smaller than 
the diffusion constant for bulk water, qualitativel y co n- 
sistent with high-T simulations of water in CNTs 17]. 

To quantitatively relate the anomalous behavior of D\\ 
with the microscopic arrangement of water molecules, we 
define ps [i] as the probability for a molecule to partici- 
pate in (exactly) i HBs, with i = 0, ... ,4, and p_p[i] as 
the probability for a molecule to have (exactly) i free 




Pressure P (GPa) 

FIG. 1: (Color online). Diffusion coefficient Dm from MC sim- 
ulations (symbols) as a function of pressure along isotherms 
for (from top to bottom) T between 972 and 483 K. For 
T < 972 K, D,| has maxima Djp x at P max (dotted-dashed 
line) and minima Djj" 111 at p mm (dashed line). Solid lines are 
from Eq. ((6| as described in the text. 



n.n. cells available for diffusion (Fig. [5]). By definition 
YjiPs[i] = J2iPF[i] = 1 and p d = 1 -ps[4] is the proba- 
bility to have defects in the HB network. Because ps [4] is 
a monotonic function also for temperatures at which D\\ 
is non-monotonic, our result show clearly that the net- 
work defects cannot be solely responsible for anomalous 
water diffusion. 

In particular, the probabilities pb[i\ behave regularly 
for any i (Fig. [5]). The probability ps[0] of a molecule 
not participating in any HB increases by increasing P, 
because of the enthalpic cost of forming HBs at high P. 
For the same reason, the other pb[i\ to form i > HBs 
decrease with P, being vanishing small for i = 3 and 4. 

Moreover, the behavior of is regular within the 

liquid phase, with the probability pp[0] to have free 
n.n. cells (high density) increasing discontinuously at the 
gas-liquid phase transition. Although pf[0] decreases by 
increasing P in the liquid phase, it is always larger than 
in the gas phase. The decrease of pf[0) is compensated 
by increases in pf[2], while pf[^\ is almost constant in 
the liquid phase and jjf[3] and pf[4] are vanishing small. 
We understand the small increase of pf[3] and Pf[4] as a 
consequence of the phase separation between liquid and 
gas, favored by the breaking of HBs due to the increase 
of P. By comparing two different temperatures, we find 
that at any P the higher T has lower pf[0], consistent 
with its lower density (Fig. [5]). 

Therefore, none of the quantities ps[i](P) or pf[i\(P) 
shows an evident non-monotonic behavior in the liquid 
phase that could be unambiguously related to the non- 
monotonic behavior of Dm(P). In the liquid phase, we 
also find a regular behavior for the average number of 
free n.n. cells around a molecule, (rip) = X)/=o / Pf\J\ 
and for the average number of HBs formed by a molecule 
(n H B> = Y,t=o b PB[ b ] (inset in Fig.©. 

We observe (Fig. [3]) that at 693 K, for P decreasing 
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FIG. 3: (Color online). Upper inset: Average numbers of 
free n.n. cells around a molecule (rip) (circles) and of HBs 
formed by a molecule (tihb) (squares) as a function of P at 
T — 693 K. The discontinuity in (np) at low P corresponds to 
the gas-liquid first-order phase transition. At higher P, both 
quantities are monotonic. Main panel: The plot of Du versus 
Wv,p along the six isotherms represented in Fig. [T] shows that 
Du is a linear function of Wv,p within the numerical noise. 
The parameters v and fj, change with T, as described in the 
text. Lower inset: Example of a cooperative (shaded) region 
of about 1 nm size, with v = 12 molecules with n.n. free cells 
(not all represented in the scheme) and (j, Vb = 5 HBs. 



FIG. 2: (Color online). Microscopic configuration changes 
along isotherms for temperatures 1006 K (squares) and 693 K 
(circles), corresponding to the onset of anomalous Du and to 
the fully developed anomaly of Du , respectively. Top: Prob- 
ability pb[i\ for a molecule to participate in i = 0, . . . , 4 HBs. 
Bottom: Probability Pf[i\ for a molecule to have i = 0, . . . , 4 
free n.n. cells available for diffusion. The discontinuity in 
Pf[i] corresponds to the liquid-gas first-order phase transi- 
tion. Among different panels, probability values change up to 
two orders of magnitude. 



from 0.4 GPa to approximately atmospheric pressure, 
the average number of HBs per molecule (hhb) increases 
from 0.3 to f=a 0.8, while the average number of free n.n. 
cells, (uf), around each molecule has a small decrease 
from 1 to 0.9. Hence, by decreasing P the change in HBs 
is large, but not the reduction of free volume around each 
molecule. 

Nevertheless, diffusion in a liquid is a process that in- 
volves more than one molecule at a time, as a conse- 
quence of the cooperative displacement of molecules, es- 
pecially at high densities [26|]. Therefore, what is rele- 
vant is not just the available volume around a molecule, 
but the volume available for diffusion Vda (dynamically 
available volume) over a region large enough to allow for 
cooperative movement of the molecules in it [261 ] . 

To find an analytic expression for quantifying the bal- 
ance of Vda an d broken HBs necessary for diffusion, we 
observe that for a microscopic configuration with aver- 



ages (n F ) and (tihb) (both functions of P and T), the 
enthalpy per molecule Eq. ([1]) is, in mean-field approxi- 
mation, 

H{P,T) = H{(n F ),{nsB)} 

= -e(4-(n F ))~{J-Pv HB )(n HB ), (3) 

where we assume that the average contribution of the 
cooperative term in Eq. ([T]) is negligible in the considered 
range of P-T, as suggested by our calculations [l9j . 

Defining Vf = ("f)/4 as the probability for each cell 
to have a n.n. cell that is free, and Vb = ("hb)/4 as 
the probability for each possible HB to be formed, the 
quantity 

W(P,T) = P F V b ^exp[-H{P,T)/(k B T)] (4) 

is the joint probability that, at a given P and T, a 
molecule forms a HB and has a free n.n. cell, with 

4 4-/ 

Z=J2PF[f]J2PBi b ] eM-H{f,b}/(k B T)} (5) 

where H{f, b} is given by Eq. ([3]). The generalization of 
W{P,T) is 

W^(P.T) = V F M 7>&~exp[-£r(P,T)/(A B T)], (6) 

and represents the probability, at given P and T, of find- 
ing v molecules with a n.n. free cell available for diffusion 



4 



within a region with /i Vb HBs (lower inset in Fig. [3]). If 
the diffusion behavior is dominated by the cooperative 
rearrangement of molecules within a given region of the 
system, then D\\ must be directly proportional to W v ^ 
for some values of v and fi. 

A comparison with our MC simulations confirms this 
hypothesis (Fig. [3]). We find that, within the numeri- 
cal noise, Dm = DW V ^ + D , where D = lA 2 /ps. The 
parameters v, fi, and Dq depend only on T, with the 
first two describing the cooperative rearranging region 
and the last associated with the case fi = 0, i. e. with 
the diffusion constant at high P. In particular, we find 
Do < 0.03 A 2 /ps, consistent with the observation that 
at very high P water recovers a normal behavior and the 
diffusion is extremely small. For the range of P-T where 
we observe the diffusion anomaly, we find v = 12.5 ± 0.5 
and n decreasing from 15 ± 1 at T = 483 K to ~ 4 ± 3 
for T > 623 K. The resulting v — 12.5 suggests that 
for the water monolayer between hydrophobic walls the 
diffusion mechanism, in the studied P-T range, requires 
a cooperative rearranging region that extends over ^3.5 
water molecules, corresponding to ~1 nm (lower inset 
in Fig. [3]). The outcome for fi conveys that at lower 
T almost all the HBs, ~18, within this cooperative re- 
arranging region need to be broken to allow for macro- 
scopic diffusion, while at higher T the number of broken 
HBs necessary for diffusion rapidly decreases, as a con- 
sequence of the reduced number of formed HBs (Fig. [2]). 

In conclusion, the result of our analysis is twofold. On 
the one hand, it clarifies the mechanisms inducing the dif- 
fusion anomaly. In particular, it shows that at constant 
T, by increasing P, the number of HBs (jihb) decreases, 
implying a decrease of the energy cost for a molecule 



to move and an increase of free volume, proportional 
to (rip), available for diffusion and competing with the 
decrease of free volume due to the increase of density. 
These competing mechanisms cause the increase of dif- 
fusion under pressurization at low P and their combined 
effect reaches a maximum at a pressure above which HB 
formation is unfavorable for enthalpic reasons. For the 
conditions considered here (high T and strong hydropho- 
bic confinement), Dy is reduced with respect to bulk. 
On the other hand, our calculations give a quantitative 
description of the phenomenon and show that the com- 
petition between free volume, available for diffusion, and 
the formation and breaking of HBs occurs within a co- 
operative region of approximately three water molecules, 
or 1 nm, in the P-T range relevant for diffusion anomaly 
in a hydrophobically nanoconfined monolayer of water. 
This result recalls the finding of Ref. [27| , where the size 
of cooperatively rearranging regions is estimated to be 
of the order 1 nm for a large number of glass forming 
materials close to the glass temperature. Here, the 1-nm 
size is found at T higher then the bulk glass tempera- 
ture, consistent with the fact that D» in confinement is 
reduced with respect to bulk. Finally, the fact that sim- 
ulations for water confined between graphite-like plates 
[TH or CNTs [H[ show a large decrease of diffusion when 
the confining length-scale is below 1 nm, suggests that 
strong confinements of water interfere with the cooper- 
ative rearrangement necessary for diffusion. Here the 1- 
nm scale emerges as a natural length associated with the 
diffusion mechanism. 
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